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ABSTRACT 


The  solution  is  obtained  for  the  reflection  and  transmission  of  time 
harmonic  plane  dilatational  or  shear  waves  at  a plane  interface  between  two 
linearly  viscoelastic  materials.  Except  in  special  cases,  the  reflected  and 
transmitted  waves  are  general  plane  waves,  i.e.,  plane  waves  whose  amplitudes 
vary  across  their  wave  fronts.  The  angles  of  reflection  and  transmission 
depend  on  the  Incident  angle  in  a more  complicated  way  than  in  the  limiting 
elastic  case  because  the  wave  speeds  of  the  reflected  and  transmitted  waves 
are,  in  general,  functions  of  the  Incident  angle.  Necessary  and  sufficient 
conditions  for  the  existence  of  Interface  waves  are  obtained.  It  is  found 
that  no  interface  waves  can  exist  for  some  materials.  For  another  class 
of  materials,  interface  waves  exist  for  discrete  angles  of  incidence.  For 
still  other  materials,  interface  waves  exist  for  all  Incident  angles  greater 

than  some  critical  angle. 
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INTRODUCTION 


This  paper  treats  the  reflection  and  transmission  of  a time  harmonic 
plane  dllatatlonal  or  shear  wave  by  a plane  interface,  y ■ 0,  between  two 
homogeneous  and  isotropic  linearly  viscoelastic  materials  as  indicated  by 
Figure  1.  The  method  of  solution  is  the  same  as  the  procedure  used  in 
solving  the  elastic  problem  by  use  of  potential  functions  except  that  it 
leads  to  the  definition  of  complex  reflection  and  transmission  angles. 

Further  analysis  of  the  formal  solution  to  obtain  the  real  reflection  and 
transmission  angles  leads  to  the  definition  of  general  plane  waves,  i.e., 
plane  waves  whose  amplitudes  vary  across  their  wave  fronts.  This  result 
is  consistent  with  previous  work*  dealing  with  the  reflection  of  plane  waves 
from  plane  rigid  and  free  boundaries  where  such  waves  were  defined  initially 
so  that  no  complex  angles  would  enter  into  the  analysis.  That  procedure 
could  also  be  applied  to  the  interface  problem;  however,  the  method  given 
here  leads  more  directly  to  a formal  representation  of  the  solution. 

The  problems  treated  here  and  in  Reference  1 have  been  considered  by 
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Lockett  by  using  a method  which  is  different  than  either  of  the  approaches 
mentioned  above.  Our  results  agree  in  general  with  his,  except  that  we 
present  a bit  more  detailed  analysis  of  the  reflection  and  refraction  phenomena 
and  obtain  a few  more  explicit  results.  In  particular,  we  obtain  necessary 
and  sufficient  conditions  for  the  existence  of  interface  waves  and  find  that, 
in  general,  such  waves  can  exist  only  for  distinct  incident  angles,  a result 
which  apparently  was  not  recognized  from  his  analysis.  We  also  obtain 
explicit  relationships  for  the  transmitted  and  reflected  wave  speeds  and 
attenuation  coefficients. 
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As  might  be  expected  from  the  previous  work,  the  phase  velocities  of 
the  reflected  and  transmitted  waves  are  functions  of  the  incident  angle 
as  well  as  material  properties,  and  the  real  reiitction  and  transmission 
angles  depend  on  the  incident  angle  in  a more  complicated  way  than  in 
the  elastic  case.  Also,  as  indicated  above,  interface  waves  are  generally 
possible  only  for  discrete  incident  angles;  whereas,  in  the  elastic  case 
or  in  a special  viscoelastic  case,  they  are  possible  for  a range  of 
incident  angles.  In  the  elastic  case  this  phenomena  is  referred  to  as 
total  reflection  or  refraction^.  In  general,  a phase  shift  between 
the  incident  and  the  reflected  or  refracted  waves  is  Introduced  by  the 
Interface . 

In  Section  I,  the  problem  is  formulated  and,  in  Section  II  the 
formal  solution,  for  either  an  incident  dilatational  or  shear  wave,  is 
obtained.  An  analysis  and  Interpretation  of  the  physical  consequences 
of  the  formal  solution  is  presented  in  Section  III.  Finally  the  results 
are  summarized  in  Section  IV. 
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I.  FORMULATION 


! 

i 


I 

[ i 


: 


% 


Let  u°f£,  t)  - u(x)  exp  {-iuit}  and  t Cx , t)  - t^(x)  «XP  denote 

the  time  harmonic  particle  displacement  vector  and  stress  tensor  where  4 > 0 
Is  the  frequency,  t Is  time,  and  s ■ (Xj,  x2,  x3)  Is  a position  vector  In 
Cartesian  coordinates.  Here,  the  subscripts  1,  j take  on  the  values  1,  2,  3. 

Hereafter,  the  complex  quantities  y,  and  will  be  referred  to  as  the  particle 
displacement  and  stress  tensor.  The  equations  of  motion,  with  body  forces 

neglected,  for  small  amplitude  time  harmonic  waves  in  a homogeneous  and  isotropic 
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linearly  viscoelastic  material  are 

( 

yAu  + (a  + y)  V (V  • u)  + pw2u  ■ 0 (1) 

where  X(u)  and  y(u>)  are  specified  complex  valued  material  properties  which 
reduce  to  the  Lam€  constants  in  the  limiting  elastic  case.  The  constant 
density  is  denoted  by  p,  and  the  gradient  and  Laplace  operators  are  denoted 
by  7 and  A.  The  constitutive  relations  are 


tjj  * + y(3u1/3Xj  + Su^/dx^ 


(2) 


where  6^  is  the  Kronecker  delta  and  1,  J ■ 1,  2,  3. 

We  consider  two-dimensional  motions  of  the  form  y * u(x,  y,  0)  where 
x^  ■ x,  *2  • y,  x}  ■ 1 are  the  usual  Cartesian  coordinates.  Since  u la 
Independent  of  z,  u3  Is  uncoupled  from  u;  and  u2  in  Eq.  1 and  u3  satlafics 
the  z component  of  Eq.  1.  The  solutions  of  this  equation  describe  shear 
motions  transverse  to  the  x,  y plane.  They  have  been  previously  Investigated 
for  a restricted  class  of  viscoelastic  materials^  and  are  not  considered  here. 


v-  * 
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Only  the  "In  plane"  notions  described  by  the  x and  y components  of  Eq.  1 are 
studied.  Hence,  we  take  u3  - 0,  Uj  - u and  u2  • v in  Eqs.  1 and  2. 

In  the  following  analysis,  the  subscript  m ■ 1,  2 is  used  to  Indicate  whicn 
medium  the  subscripted  variable  refers  to,  i.e.,  m » 1 implies  y > 0 and  m - 2 
implies  y < 0,  as  shown  in  Fig.  1.  The  subscript  n • 1,  2 is  used  to  denote 
dilatatlonal  and  shear  waves  respectively,  and  the  subscript  t • 1,  2 denotes 
the  case  where  a dilatatlonal  (t  ■ 1)  or  a shear  (l  • 2)  wave  is  Incident  at  the 
Interface. 

Eight  potential  functions,  ♦ *t*  defined  such  that,  for  l,  m * 1,  2, 


ut»  “ 3*tmi/3x  + 3*lm2/3y’ 


3*  /3y  - 3*. 


(3) 


’ tm  -'1ml' tm2' 

are  the  displacement  components.  The  displacement  vector  obtained  from  Eq.  3 
is  a solution  of  Eq.  1 provided 


+ k2  ♦.  - 0, 

imn  mn  tmn 


1,  m,  n a 1,  2, 


(4) 


where 


k«, ! “/s™; 


•l 


<*«  + 2wm>/p»;  Si2  " Mm/Pm* 


(5) 


Here,  k is  the  "complex  wave  number"  for  dilatatlonal  (n  **  1)  or  shear  (n  “ 2) 

mn 

waves  in  medium  m which  may  be  expressed  as1 


(6) 


for  m,  n - 1,  2.  The  complex  quantities,  , are  referred  to  a9  "complex 

wave  speeds"  because  of  their  analogy  to  the  elastic  wave  speeds.  It  will  be 

shown  that  the  real  quantities,  C , are  least  upper  bounds  for  the  dllatational 

mn 

(n  - 1)  and  shear  (n  ■ 2)  wave  speeds  In  medium  m.  The  requirement  that 

0 < ft  < h/2  comes  from  the  condition  that  plane  waves  do  not  grow  in  amplitude 
ran 

as  they  propagate1. 

An  effort  has  been  made  to  keep  the  notation  consistant  with  that  in 
Ref.  1.  However,  in  order  to  avoid  confusion  with  the  Kronecker  delta,  0 is 
used  for  the  variable  defined  as  6 in  the  previous  work.  Also,  we  have 
reseived  $ to  denote  an  angle.  Hence,  the  potential  functions  are  here  denoted 
by  ♦. 

The  normal  stresses  in  the  x and  y directions,  for  medium  m and  incident 

wa*e  i,  are  denoted  by  a . and  o . , and  the  corresponding  shear  stress  is 

x x>ra  y im 

denoted  by  x.  . Thus,  from  Eqs.  2 and  3,  for  l,  m ■ 1,  2, 


°xtm  “ 2‘i«(5,2*t.l/3x2 


/ i>xdy) 


V.  ‘ 2".<,!V/a,a  ' - x.  “ii  *t.  • 
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The  boundary  conditions  at  the  interface,  y ■ 0,  are,  for  i ■ 1,  2, 


I 


utl(x,  o)  - u^x,  0)1 

vn(x,  0)  ■<  V 

y*l(X’  0)  “ V2(X*  0)’ 

Ttl(x*  °)  “ T 

(9) 
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The  total  potential  functions  are  defined,  for  l,  n • 1,  2,  by 


*lln  “ 6tn  % + *Un;  *t2n  “ *l2n 


(10) 


where  6.  is  the  Kronecker  delta.  The  incident  wave,  Y.,  is  a dilatatlonal 

Jin  * 

wave  if  l ■ 1 and  a shear  wave  if  1*2.  The  reflected  and  refracted  waves 

are  denoted  by  4,  . The  potential  functions  for  l,  a,  n ■ 1,  2 are  assumed 
imn 

to  be  of  the  form 


Ti  " exp  * 

iji.  ■ (l.R.  S /w)  exp  (ik  (r„  *x))  , 

Imn  l imn  mn  mn  fcmn  3 


(11) 


where 


lt  2 i sin  et  - £ cos  0£,  x = ix  + 1y, 


(12) 


- tmn  5 i 8ln  Ct»n  + J7m  008  W Em  5 (_1) 


m + l 


Here,  i,  j are  unit  vectors  along  the  positive  x and  y axes.  The  angle,  6^,  is 
the  angle  between  the  incident  dilatatlonal  (t  ■ 1)  or  shear  (i  “ 2)  ray  r^ 
and  j.  The  angles  ( are,  in  general,  complex.  If  they  are  real,  then  they 
are  between  f and  J for  m ■ 1 and  between  and  *or  ® ■ 2 (and  are  given 

b*  Ctmn  “ 0tmn) * S*8  Flg‘  U 


- ;>v  ^7  * 

1'*—  W * 

■ * * t 

A **  ,**  t 


Ic  may  be  shown  that  1^  is  the  displacement  amplitude  of  the  Incident 

wave,  normal  to  Its  wavefront,  at  the  origin  x - y » 0.  Similarly,  R.  I,  is 

x.mn  Jt 

the  displacement  amplitude  of  the  reflected  or  refracted  wave,  normal  to  its 
wavefront,  at  the  origin.  Thus,  , m,  n ■ 1,  2,  are  the  displacement  (or 
velocity)  reflection  and  refraction  coefficients  for  1 ■ 1,  2. 

Note  that  the  incident  wave  is  a plane  wave  in  the  usual  sense,  l.e.,  it 
has  constant  amplitude  across  its  wavefront.  When  C.  is  complex,  it  will  be 
shown  that  the  reflected  or  refracted  wave  is  a "general  plane  wave",  l.e.,  it 
is  attenuated  across  its  wavefront. 

From  Eqs . 3,  4,  11  and  12  it  may  be  shown  that 


utm 

ikml  8in  <lml 

ikm2  em  Coa  51«2 

vim 

e»lk»l  co8  Ctml 

“lkm2  8ln  «im2 

°xlm 

- 

-k2 . (A  + 2u  sin2 
ml  m m 

'■tml) 

“umkm2em  8ln  2<im2 

Q 

*< 

a 

-k2,(X  + 2u  cos2 

ml  m m 

l’tml) 

umk»2em  8ln  2i*tm2 

Tlm 

-e  u k2  sin  2 4, 
m m mi  tmi 

~umkm2  cos  2Cim2 

■*  MB. 

— 

+ 6ml 


i(6tlkU  810  91  - 6i2k12  COS  62> 


-i(6tlkll  co»  ®1  + 6t2k12  8ln  ®2> 
Sil 


"6llkll(Xl  + 2wl  8ln2  61>  + 6t2Wlk12  8ln  262 


<X1  + 2m1  COs2  01)  * aC2Ulk12  8ln  202 

6UWlk?i  810  26l  - ®t2**lk12  608  262 
8 


1ml 


lm2 


(13) 
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II.  SOLUTION 


The  boundary  conditions  given  by  Eq.  9 require  that  all  of  the 
exponentials  in  Eq.  13  be  equal  on  y ■ 0.  Hence,  we  obtain  the  complex 
Snell's  law,  for  1,  m,  n * 1,  2, 


sin  c. 


tarn 


S 

mn 


8ln  Vs  u • 


(14) 


where  the  complex  wave  speeds,  S , are  defined  by  Eq.  5.  Rearranging  the 

mn 

expressions  resulting  from  the  boundary  conditions  and  combining  the  result 
with  Eq.  14,  we  obtain  the  linear  set  of  equations 


(15) 


where 


sin  C 


cos  c 


ill 


ill 


?1  * 


11 


-sln0 


CO£0. 


PjSjjCos  2; 


cos  c 
-sin  c 


°1  SllCOS  <112 

PlS 

°1  Si28ln24iu 

PtS 

1 12 


p.S2  sin  20 
I 12  1 


112 

-8iD  <121 

cos  c. 

112 

C°8  <121 

sin  c 

"ln  2<112 

P2S21cos 

2<122 

°2S22 

C°8  2C112 

P2S228ln 

S21 

2<121 

_P2S22 

r 

r 

1 

COS0. 


sine. 


122 

122 


'122 


’*22 


-PjSjjSin  202 


-PiS12cos  202 


• *1 


11 1 


112 


121 


122 


(17) 
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If  det  1^1  i 0,  then  the  reflection  and  refraction  coefficients, 

R.  , are  determined  by  Inverting  Eq.  15.  In  general,  R is  complex. 

£ntn  jcinn 

Hence,  there  is  a phase  shift  between  the  Incident  wave  and  the  reflected 

and  refracted  waves  in  general.  In  order  for  the  formal  solution  obtained 

above  to  be  unique,  unique  values  for  i,ra,n  ■ 1,2,  must  be  specified 

(since  the  inversion  of  Eq.  14  yields  multivalued  roots  for  the  complex 

angles,  C,  ).  In  the  following  section,  the  physics  of  the  reflection  and 
£mn 

reflection  phenomena  will  be  studied,  and  the  information  presented  there 
will  allow  a unique  inversion  of  Eq.  14. 
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III.  PHYSICAL  INTERPRETATION  OF  THE  SOLUTION 
Let  the  complex  angle,  5 n>  be  represented  by 

^imn  a£ran  + *^imn 


(18) 


where  a.  and  8„  real  numbers.  Combining  Eqs.  5,  6,  1A,  and  18;  noting 
Jcmn  x.mn 

that  sin  - sin  a„  cosh  8„  + i cos  a.  sinh  6.  ; and  equating 

Jtmn  Jtmn  Jtnin  x.mn  imn 

the  real  and  imaginary  parts  of  the  resulting  expression,  we  obtain 


sin  a.  cosh  8„  - ^ . cos  A„ 

imn  imn  imn  imn  , 


cos  a.  sinh  6.  ■ T sin  A. 

imn  imn  imn  imn 


(19) 

(20) 


where 


r.  = y cos  0 / cos  Olo  > 0, 

imn  imn  mn  u - 


Yimn  S Cmn  8ln  ®i  /Cli  - 0 ’ 


^imn  ^lj.  ^mn 


(21) 

(22) 

(23) 


Combining  Eqs.  19  and  20  and  using  simple  identities,  we  obtain 


sinh2 

6.  + sin2 

imn 

°imn 

r 

imn  , 

(2A) 

2 

cosh2 

8.  - cos2 

Jcmn 

a m 

imn 

r 

imn  , 

(25) 

8ln  “tnm 

Jtmn 

cosh 8.  sin 
Linn 

^imn 

cos  a.  sinh  8.  cos  A. 

imn  imn  imn 

(26) 

11 


ft 


Combining  Eq.  24  with  the  square  of  Eq.  26  and  noting  that 


cosh28„  - sinh28,  ■ 1,  we  obtain  a quadratic  in  sinh2B.  whose 
imn  Iran  tmn 

solution  for  sinh28„  > 0 is 
Jcmn  ~ 


sinh20lmn  = 2 {ri 


- 1 + 


[2  , 2 

(1  - T )2  + 4T  sin2A. 
imn  tmn  1 


Note  that  4.  is  real  if  and  only  if  6.  m 0 in  Eq.  18.  Thus,  from 
Jcmn  Jcmn 

Eq.  27,  necessary  and  sufficient  conditions  for  to  be  a real,  non- 
zero angle  are  A.  =0  and  0 < T.  <1.  In  this  case,  4.  * a.  * 9.  is 

imn fcmn  ~ imn  imn  imn 

given  by  Eq.  14  or  Eq.  19.  Note  that  rimn  = 0 (0^=  0)  is  a sufficient 

condition  for  8.  * 0,  and  in  this  case  a.  ■ 0 from  Eq.  19.  Hence, 

Jcmn  Jcmn 

if  the  incident  wave  is  normal  to  the  interface,  all  reflected  and 

refracted  waves  propagate  in  directions  normal  to  the  interface. 

If  8.  ^ 0,  then  the  proper  sign  for  sinh  6.  * Isinh  8,  I given 

Jcmn  Jcmn  Jcmn 

by  Eq.  27  must  be  determined  before  the  solution  is  complete.  Since 

r 2 0,  Eq.  19  indicates  that  (note  that  0 < fl  < it/2,  so-ir/2  < A.  < «/2) 
x,mn  mn  Jcmn 


0 < O.  < Tf 
- {.mn 


Eq.  20  then  indicates  only  that  cos  a.  and  8.  either  have  the  same 

Jcmn  Jcmn 

or  the  opposite  sign  depending  on  the  sign  of  ^imn*  Thus,  the  form  of 

the  information  presented  so  far  is  not  sufficient  to  determine  the  proper 

sign  for  8„  . We  shall  obtain  the  necessary  information  by  considering 

JCmn 

the  real  properties  of  the  reflection  phenomena. 


Combining  Eqs.  11,  12,  and  18  and  simplifying,  we  obtain 


ikmn^mn  * (lw/cimn)  ( + 1 D*mnW 


where 


c.  = C 
Ian  mn  Imn  , 


= Tl  + sinh2  8 sec2fl  "I 
imn  1 imn  mn  I 


p,  = i sin  0.  + je  cos  8.  , 

-tmn  - Imn  - m tmn 


?imn  = i sin  (almn  + nimn)  + j-emCos(a*mn  + nimn)  * 


0imn  “ almn  ^Imn  , 


tan  *imn  “ tan  Qmn  tanh  Bimn  * ~1'/2  - **mn  - w/2  * 


tan  n - cot  8 tanh  6.  , -it/2  < n „<  */2, 

JliTjTi  mil  i.  mil  • mn 


°imn  ” ^imn 


[tan2  8 

' 


cosh28._  + sinh2  g.  1 ** 


mn  imn 


In  Eq,  29,  c is  the  real  wave  speed,  p and  a.  are  real  unit 
tmn  __  tmn  - tmn 

vectors* in  the  directions  of  propagation  and  damping,  and  is  the 

total  damping  coefficient.  From  Eq.  30,  £.  is  the  ratio  of  the  wave 

xnin 

speed  to  the  material  property  C given  by  Eq.  7,  and  from  Eq.  31 

mn 

0 < C.  < 1.  Hence,  the  material  properties  C are  least  upper  bounds 

xmn  . mn 
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on  the  wave  speeds.  Note,  from  Eq.  31,  that  £ 


1 If  and  only 


tmn 

if  B.  ■ 0.  Thus,  a necessary  and  sufficient  condition  for  waves 
imn  1 

to  propagate  with  speeds  specified  entirely  by  material  properties. 

i.e.  independent  of  the  incident  angle,  is  for  to  be  real.  Note 

that  this  is  always  the  case  at  normal  incidence  since  - 0 in  that 

tmn 

case. 

The  propagation  and  damping  unit  vectors,  p.  and  a.  are  represented 

- xmn  - Jcmn 

schematically  in  Figure  2 for  m - 1.  Note  that  they  coincide  if  and 

only  if  4> „ - n«  which  occurs  if  and  only  if  8.  ■ 0.  Thus, 

inn  Junn  JCmn 

a necessary  and  sufficient  condition  for  the  surfaces  of  constant  phase. 

p.  • x ■ constant,  and  surfaces  of  constant  amplitude,  a,  . x ■ constant, 

- xron  - mn  - 

to  be  parallel  is  for  c.  to  be  real.  Otherwise,  the  waves  are  attenuated 

*nm 

in  directions  other  than  the  direction  of  propagation. 

The  real  reflection  or  refraction  angle,  6.  , will  now  be 

Jcmn 

determined.  Combining  Eqs.  34,  35,  and  31  to  obtain  an  expression  for 


sin  » and  then  comparing  the  result  with  the  sum  of  Eq.  19  and 

tan  U times  Eq.  23,  we  find 
mn 


9in  9 tmn  ^imn  Ytmn  ~ 


(38) 


which  explicitly  determines  the  reflection  and  refraction 


This  result 


is  found  from  Eqs.  31  and  27,  is  precisely  the  s< 


angles  once  ( 


form  obtained  for  the  reflection  of  plane  waves  from  free  and  rigid 


y <1  since  the  analysis  has  required 


boundaries 


that  0.  be  real 
tmn 


Cl 

* 


The  attenuation  vector,  a.  , may  be  decomposed  Into  components 

*■  Xmn 

along  and  normal  to  Plmn  by  noting  that  a^  - (a^  • p^)  Ptmn  + 

(a„  • t.  ) t.  where  t„  * k x p.  and  k is  a unit  vector  in  the 

-imn  - imn  ~ tmn  -imn  . - torn 

- z direction.  Hence,  the  unit  vector  tangent  to  the  reflected 
or  refracted  wave  front,  positive  when  pointing  away  from  the  Interface, 
is 


t.  - -i  cos  0 + je  sin  0.  . 

- xmn  - xmn  m xmn 


Hence,  Eq.  29  may  be  transformed  to 


(39) 


lkInr|  ( r • *)  “ (i»/c  ) 

mn  - xmn  - xmn 


h 


1 CLn  tanQmn)  (?tmn  * 5> 


-id_  sec2  n sinh  Bcooh  6 (t 
xmn  mn  xmn  xmn  - xmn 


. x)*J 


(40) 


Hence,  the  longltudional  damping  coefficient  (in  the  direction 


°f  ?£mn)  18 


dimn  1 (w/Cmn}  *imn  tM  “mn  * 


(41) 


and  the  transverse  damping  coefficient  (in  the  direction  of  ) is 


T.  2 - (w/C  ) C.__  sec2  a sinh  B.  cosh  B 
xmn  mn  xmn  mn  xmn 


4mn 


(42) 


Thus,  the  wave  is  transversely  damped  away  from  the  interface  if  3tmn<  0 

and  toward  the  interface  if  B„  > 0.  Note  also  that  T..  • 0 if  and 

xmn  xmn 

only  if  Btmn  * 0.  Thus,  the  significance  of  the  sign  of  sinh  0Jmn  now 
becomes  clearer. 
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Multiplying  Eq.  19  by  tan  adding  the  result  to  Eq.  20, 

simplifying  and  noting  Eqs.  34,  35,  38,  31  and  42,  we  obtain 


<i>r 


imn 


^ ^imn 


tan  0 


mn 


tan  (1 


U 


] - c. 


mn  imn 


cos  6 


tmn 


(43) 


We  note  that  Eq.  43  is  consistent  with  the  previous  findings  on  necessary 

and  sufficient  conditions  for  8.  *0  (and,  therefore,  T.  ■ 0)  in  that 

x,mn  Jcmn 

those  conditions  are  Included  in  the  set  of  conditions  for  which  the  left 

side  of  Eq.  43  vanishes.  From  the  physical  situation,  |6tmnl  £ */2.  Thus 

Eq  . 28  indicates  that  0 < 0.  < rr/2  so  that  cos  6.  is  non-negative. 

^ imn  - imn 

Hence,  the  signsof  T.  and  8.  i 0 are  determined  by  the  left  side  of 

xmn  xmn 

Eq . 43  provided  both  sides  are  not  identically  zero.  Excluding,  for  the 

moment,  interface  waves  (6.  - n/2)  and  the  case  of  no  transverse  damping, 

x,mn 


8lgn  Timn  " 8lgn  (~Bimn)  “ 8ign  ULn  tan  °mn  " tan  flU)! 


if 


imn 


i o, 


0imn  * "/2* 


(44) 


Note  that  if  ft  > ft  , then  T.  <0  since  <1  and  the  reflected 

1 1,  mn  tmn  tmn  — 

or  transmitted  wave  is  transversely  damped  toward  the  interface. 

We  now  determine  necessary  and  sufficient  conditions  for  Interface 

waves  to  exist,  l.e.,  6.  - tt/2.  To  facilitate  the  argument,  we  consider 

inn 

three  cases:  a)  0 j*  £},  , b)  £1  ■ £1  f 0,  and  c)  £1  ■ £1  « 0. 

mn  It  mn  lj,  mn 

If  £1  + (1  , then  from  Eq.  28  8.  + 0.  Furthermore,  we  have 

mn  1 £ Jtmn 

already  shown  only  normally  reflected  and  transmitted  waves  exist 
if  • 0 (normal  incidence).  Hence,  from  Eq.  43,  a necessary  and 
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sufficient  condition  for  6 


fmn 


n/2  (If  ft  t ft  ) is 
inn  „ i 


(,*'  tan  ft  * tan  ft 
x.mn  mn  1 ^ 


(4ii 


Note  that  this  condition  can  not  be  satisfied  11  eithei  ft  » 0 or 

mn 

ft  *■  0,  or  If  Q » ft  (since  i,  < 1),  Thus,  no  interface  wave  can  exist 
1 1_  1 1 mn 

In  these  cases-  Hence,  if  ft  + 0 and  ft.,  t*  0 and  li  ft  i ft,.,  then  from 

mn  » it  mn  i Jt 

Eg.  38.  an  interface  wave  exists  if  and  only  If 


t(  = tan  ft  /tan  ft.  < 
Vmn  mn  1 i - 


i,  (i  * ft. 

mn  l 


l ’ 


(46) 


Note  •’hat  li  an  interface  wave  occurs,  then  it  occurs  for  a discrete  incident 
angle 

If  ft  » ft.  . f 0,  then  A.  “0  and  we  have  previously  found  that 
mn  1 k.  imn 


6.  * 0 if  and  only  if  T,  ■ y.  < 1.  (See  Eq.  28.)  Thus,  from  Eq.  31, 

^f.mn  “ 1 in  this  case.  Hence,  from  Eq.  38,  an  interface  wave  is  possible 

if  and  only  if  Yimn  > 1.  In  fact,  from  Eqs.  43  and  38,  a necessary  and 

sutiicient  condition  for  9.  ■ »/2  (if  ft  - ft,  + 0)  is 

Imn  mn  1 1 


Yimn  " l' 


(47) 


To  make  this  clear,  note  that  the  left  side  of  Eq.  43,  in  this  case  reduces 

to  r ( 1- ) tan  ft  , which  can  vanish  if  and  only  if  - 1.  Thus, 

iron  imn  mn  imn 

from  Eq , 38,  9.  ■ n /?  if  and  only  if  y.  “ 1.  Note  that  this  implies 

imn  Imn 

that  an  interface  wave  can  exist  for  only  a distinct  incident  angle. 

If  ft  - ft,.  - 0,  then  8.  i 0 if  v.  > 1 as  in  the  previous  case, 
mn  1 JL  imn  imn 

Further,  the  left  side  of  Eq.  43  is  identically  zero.  Thus,  if  8^  / 0, 


' V*  _*  - 


17 


then  0.  ■ w/2.  Note  that  £«  *1  for  y.  < 1.  Thus,  from  Eq.  38, 

f ran  v mn  v mn  — 


"£mn 


imn- 


0£mn  < 1,72  for  Y£mn  < and  0£mn  " 11/2  if  Y£mn  " K HenCC’  0£mn  “ "/2 

for  y«  > 1.  For  the  reciprocal  argument,  note  that  if  6„  - n/2,  then 

x.mn  -*  v 


tmn 


from  eq  38,  y * 1.  Then  since  0 < r < 1,  it  is  clear  that 

Jcmn  £ran  Jtz&n  — 

Yimn  - a necessary  and  sufficient  condition  for  9g  ■ it/2 


imn- 


(if  n - 0,  - 0)  is 

mn  t i 


Yimn  - 


(48) 


This  implies  that  an  interface  wave  exists  for  any  incident  angle  such  that 


sin  6.  > C,  /C 

£ 1 1 mn 

It  remains  to  specify  the  proper  sign  for  T.  and  8.  when  0, 

x,mn  Jcmn  Jcmn 

In  all  cases  where  Interface  waves  exist,  we  choose 


£nm 


> 0,  6.  <0  for 

Jcmn  — 


%«,  ■ 'n • 


(«) 


so  that  the  solution  Is  bounded  for  increasing  y.  This  choice  also  provides 

consistency  with  the  limiting  elastic  case.  Note  that  if  0 - ft, » i 0, 

mn  * *’ 

then  from  Eqs.  47  and  28,  0.  ■ 0. 

jcmn 
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IV.  DISCUSSION 


From  the  preceding  analysis,  we  can  make  several  general  observations 

about  the  reflection  and  transmission  phenomena  at  the  Interface.  In  general, 

the  reflected  and  transmitted  waves  are  general  plane  waves  (both  longitudinal 

and  transverse  damping)  whose  wave  speeds  and  attenuation  coefficients  are 

functions  of  the  incident  angle  as  well  as  material  properties.  The  least 

upper  bound  for  the  wave  speed  Is  the  material  property  C where  m,  n - 1,2. 

mn 

In  general,  interface  waves  exist  for  distinct  angles  of  Incidence  if  they 
are  possible  at  all.  In  general,  there  is  a phase  shift  between  the  incident 
and  reflected  or  transmitted  wave. 

The  exceptions  to  the  above  general  statements  center  around  special 

cases  for  the  complex  wave  speeds  S . If  ImS  - ImS,.,  then  the  reflected 

mn  mn  1 * 

or  transmitted  wave,  ♦ . , is  a plane  wave  in  the  usual  sense  (constant 

amplitude  across  the  wave  front)  unless  it  is  an  interface  wave,  in  which 

case  it  becomes  a general  plane  wave.  If  ImS^  ■ ImSj^  i 0,  an  interface 

wave  can  exist  only  for  distinct  incident  angles;  but  if  ImS  ■ ImS  ■ 0, 

then  interface  waves  exist  for  any  incident  angle  greater  than  some  critical 

angle.  Also,  the  only  case  where  the  wave  speed  is  Independent  of  the 

incident  angle  is  when  ImS  ■ ImS... 

mn  1 1 

Another  special  case  of  some  interest  li  0..  > fl . In  this  case, 

no  Interface  wave  is  possible.  Also,  It  was  shown  that  the  reflected  or 
transmitted  wave  Is  attenuated  across  its  front  toward  the  interface  in 
this  case.  This  result  is  contrary  to  the  elastic  case  of  internal  reflection 
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where  the  wave  is  attenuated  away  from  the  Interface.  For  the  reflected 

2 

waves  (m  - 1) , Lockett  has  suggested  that  when  the  incident  wave  is  a 
shear  wave  (2^  > flln  for  most  real  materials.  The  transmitted  wave 
properties  cannot  be  compared  without  specifying  specific  materials. 

Hence,  no  such  relationship  between  and  ft2[)  generally  applicable. 
The  results  presented  here  are  entirely  consistent  with  those  for 

rigid  and  free  boundaries  in  Reference  1.  In  generalities,  these  results 

2 

also  agree  with  Lockett  s results  . However,  there  is  at  least  one  point 
where  an  inconsistency  occurs.  Lockett  indicates  that  if  interface  waves 
exist,  they  occur  for  incident  angles  greater  than  some  critical  angle. 

It  was  shown  here  that,  except  for  the  above  mentioned  special  cases, 
interface  waves  occur  for  discrete  incident  angles,  i.e.,  the  reflected  or 
transmitted  ray  moves  away  from  the  Interface  if  the  Incident  angle  Is 
Increased  beyond  the  critical  angle. 

The  aolutlon  presented  here  includes  all  cases  of  homogeneous  and 
Isotropic  linearly  viscoelastic  materials  including  the  limiting  elastic 
case.  Once  the  complex  angles  C&mn  are  determined,  the  results  from 
Section  II  can  be  applied  to  determine  the  field  for  a fixed  frequency  w. 
If  pulse  propagation  is  of  interest,  then  a Fourier  synthesis  may  be  used 
to  solve  the  transient  problem. 


i 


REFERENCES 


1.  H , F.  Cooper,  Jr.,  and  E.  L,  Reiss,  J,  Acoua*  , Soc  Am  J9,  1133- 
1138  (1966). 

2 F,  J,  Lockett,  J.  Mech.  Phys.  Solids,  10,  53-64  (1962). 

3.  H.  Kolsky,  Stress  Waves  in  Solids.  (Dover  Publications,  Inc., 

New  York,  1963) . 

4.  D,  R . Bland,  Theory  of  Linear  Viscoelasticity.  (Pergamon  Press,  Inc., 
New  York , 1960) . 

5.  H T„  O'Neil,  Phys,  Rev,  75,  928  - 935  (1949). 


21 


l— - i. 


ki 


LIST  OF  FIGURES 


Figure  1 - Reflection  and  Refraction  Phenomena  at 
the  Interface  y ■ 0. 

Figure  2 - Propagation  and  Damping  Directions  of 

the  Reflected  Wave  i|>. 

tin 


22 


FIGURE  I - REFLECTION  AND  REFRACTION  PHENOMENA  AT  THE 
INTERFACE  Y»0 
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FIGURE  2-  PROPAGATION  AND  DAMPING  DIRECTIONS  OF 
REFLECTED  WAVE  ^(f) 


